Continuity equation

A continuity equation in physics is an equation that describes the transport of a conserved quantity. Since mass, energy, momentum, electric charge and other natural quantities are conserved under their respective appropriate conditions, a variety of physical phenomena may be described using continuity equations.

Continuity equations are the (stronger) local form of conservation laws. All the examples of continuity equations below express the same idea, which is roughly that: the total amount (of the conserved quantity) inside any region can only change by the amount that passes in or out of the region through the boundary. A conserved quantity cannot increase or decrease, it can only move from place to place.

Any continuity equation can be expressed in an "integral form" (in terms of a flux integral), which applies to any finite region, or in a "differential form" (in terms of the divergence operator) which applies at a point.

Contents

General equation

Preliminary description

As stated above, the idea behind the continuity equation is the flow of some property, such as mass, energy, charge, momentum, and even probability, through surfaces from one region of space to another. The surfaces, in general, may either be open or closed, real or imaginary, and have an arbitrary shape, but are fixed for the calculation (i.e. not time-varying, which is appropriate since this complicates the maths for no advantage). Let this property be represented by just one scalar variable, q, and let the volume density of this property (the amount of q per unit volume V) be φ, and the all surfaces be denoted by S. Mathematically, φ is a ratio of two infinitesimal quantities:

\varphi = \frac{{\rm d} q}{{\rm d} V},

which has the dimension [quantity][L]-3.

There are different ways to conceive the continuity equation:

  1. either the flow of particles carrying the quantity q, described by a velocity field v, which is also equivalent to a flux f of q (a vector function describing the flow per unit area per unit time of q), or
  2. in the cases a velocity field is not useful or applicable, the flux f of the quantity q only (no association with velocity).

In each of these cases, the transfer of q occurs as it passes through two surfaces, the first S1 and the second S2.

The flux f should represent some flow or transport, which has dimensions [quantity][T]−1[L]−2. In cases where particles/carriers of quantity q are moving with velocity v, such as particles of mass in a fluid or charge carriers in a conductor, f can be related to v by:

\mathbf{f} = \varphi \mathbf{v} .

This relation is only true in situations where there are particles moving and carrying q - it can't always be applied. To illustrate this: if f is electric current density (electric current per unit area) and φ is the charge density (charge per unit volume), then the velocity of the charge carriers is v. However - if f is heat flux density (heat energy per unit time per unit area), then even if we let φ be the heat energy density (heat energy per unit volume) it does not imply the "velocity of heat" is v (this makes no sense, and is not practically applicable). In the latter case only f (with φ) may be used in the continuity equation.

Elementary vector form

Consider the case when the surfaces are flat and planar cross-sections. For the case a velocity field can be applied, dimensional analysis leads to this form of the continuity equation:

 \varphi_1 \bold{v}_1 \cdot \bold{S}_1 = \varphi_2 \bold{v}_2 \cdot \bold{S}_2

where

Notice the dot products  \bold{v}_1 \cdot \bold{S}_1, \, \bold{v}_2 \cdot \bold{S}_2 \,\! are volumetric flow rates of q. The dimension of each side of the equation is [quantity][L]-3·[L][T]-1·[L]2 = [quantity][T]-1. For the more general cases, independent of whether a velocity field can be used or not, the continuity equation becomes:

 \bold{f}_1 \cdot \bold{S}_1 = \bold{f}_2 \cdot \bold{S}_2

This has exactly the same dimensions as the previous version. The relation between f and v allows us to pass back to the velocity version from this flux equation, but not always the other way round (as explained above - velocity fields are not always applicable). These results can be generalized further to curved surfaces by reducing the vector surfaces into infinitley many differential surface elements (that is S → dS), then integrating over the surface:

 \iint_{S_1} \varphi_1\bold{v}_1 \cdot {\rm d}\bold{S}_1 = \iint_{S_2} \varphi_2\bold{v}_2 \cdot {\rm d}\bold{S}_2

more generally still:

 \iint_{S_1} \bold{f}_1 \cdot {\rm d}\bold{S}_1 = \iint_{S_2} \bold{f}_2 \cdot {\rm d}\bold{S}_2

in which

N.B: the scalar area S and vector area S are related by \mathrm{d}\mathbf{S} = \mathbf{\hat{n}}\mathrm{d}S. Either notations may be used interchangeably.

Differential form

The differential form for a general continuity equation is (using the same q, φ and f as above):

\frac{\partial \varphi}{\partial t} %2B \nabla \cdot \mathbf{f} = \sigma\,

where

This general equation may be used to derive any continuity equation, ranging from as simple as the volume continuity equation to as complicated as the Navier–Stokes equations. This equation also generalizes the advection equation. Other equations in physics, such as Gauss's law of the electric field and Gauss's law for gravity, have a similar mathematical form to the continuity equation, but are not usually called by the term "continuity equation", because f in those cases does not represent the flow of a real physical quantity.

In the case that q is a conserved quantity that cannot be created or destroyed (such as energy), this translates to σ = 0, and the continuity equation is:

\frac{\partial \varphi}{\partial t} %2B \nabla \cdot \mathbf{f} = 0\,

Integral form

By the divergence theorem (see below), the continuity equation can be rewritten in an equivalent way, called the "integral form":

 \frac{{\rm d}q}{{\rm d}t} %2B \iint\limits_{S}\!\!\!\!\!\!\!\!\!\!\!\subset\!\supset (\mathbf{f}\cdot\mathbf{\hat{n}})\,{\rm d}S = \Sigma

where

In a simple example, V could be a building, and q could be the number of people in the building. The surface S would consist of the walls, doors, roof, and foundation of the building. Then the continuity equation states that the number of people in the building increases when people enter the building (an inward flux through the surface), decreases when people exit the building (an outward flux through the surface), increases when someone in the building gives birth (a "source" where σ > 0), and decreases when someone in the building dies (a "sink" where σ < 0).

Derivation and equivalence

The differential form can be derived from first principles as follows.

Derivation of the differential form

Suppose first an amount of quantity q is contained in a region of volume V, bounded by a closed surface S, as described above. This is equal to the amount in already in V, plus the generated amount s (total - not per unit time or volume):

q = \int \mathrm{d}q = \iiint_V \varphi \mathrm{d}V %2B s.

where s is such that:

\Sigma = \frac{\mathrm{d}s}{\mathrm{d}t}

The rate of change of q leaving the region is simply the time derivative:

\begin{align} \frac{\partial q}{\partial t} & = -\frac{\partial }{\partial t} \iiint_V \varphi \mathrm{d}V %2B \frac{\partial s}{\partial t} \\
& = - \iiint_V \frac{\partial \varphi}{\partial t} \mathrm{d}V %2B \Sigma ,\\
\end{align}

where the minus sign has been inserted since the amount of q is decreasing in the region. (Partial derivatives are used since they enter the integrand, which is not only a function of time, but also space due to the density nature of φ - differentiation needs only to be with respect to t). The rate of change of q crossing the boundary and leaving the region is:

\frac{\partial q}{\partial t} = \iint\limits_{S}\!\!\!\!\!\!\!\!\!\!\!\subset\!\supset (\mathbf{f}\cdot\mathbf{\hat{n}})\mathrm{d}S ,

so equating each expression for  \partial q / \partial t:

\iint\limits_{S}\!\!\!\!\!\!\!\!\!\!\!\subset\!\supset (\mathbf{f}\cdot\mathbf{\hat{n}})\mathrm{d}S = - \iiint_V \frac{\partial \varphi}{\partial t} \mathrm{d}V %2B \Sigma,

Using the divergence theorem on the left-hand side:

\iiint_V (\nabla\cdot\mathbf{f}) \mathrm{d}V = - \iiint_V \frac{\partial \varphi}{\partial t} \mathrm{d}V %2B \iiint_V \sigma \mathrm{d}V.

This is only true if the integrands are equal, which directly leads to the differential continuity equation:

\begin{align} & - \frac{\partial \varphi}{\partial t} %2B \sigma = \nabla\cdot\mathbf{f}, \\
& \nabla\cdot\mathbf{f} %2B \frac{\partial \varphi}{\partial t} = \sigma  \rightleftharpoons \nabla\cdot(\varphi \mathbf{v}) %2B \frac{\partial \varphi}{\partial t} = \sigma.\\
\end{align}

Either form may be useful and quoted, both can appear in hydrodynamics and electromagnetism, but for quantum mechanics and energy conservation, only the first may be used. Therefore the first is more general.

Equivalence between differential and integral form

Starting from the differential form which is for unit volume, multiplying throughout by the infinitesimal volume element dV and integrating over the region gives the total amounts quantities in the volume of the region (per unit time):

\begin{align} & \iiint_V \frac{\partial \varphi}{\partial t} \mathrm{d}V %2B \iiint_V (\nabla \cdot \mathbf{f}) \mathrm{d}V = \iiint_V \sigma \mathrm{d}V \\
& \frac{\mathrm{d}}{\mathrm{d} t} \iiint_V \varphi\mathrm{d}V %2B \iiint_V (\nabla \cdot \mathbf{f}) \mathrm{d}V = \Sigma \end{align} \,\!

again using the fact that V is constant in shape for the calculation, so it is independent of time and the time derivatives can be freely moved out of that integral, ordinary derivatives replace partial derivatives since the integral becomes a function of time only (the integral is evaluated over the region - so the spatial variables become removed from the final expression and t remains the only variable).

Using the divergence theorem on the left side obtains the integral form:

\frac{\mathrm{d} q}{\mathrm{d} t} %2B \iint\limits_{S}\!\!\!\!\!\!\!\!\!\!\!\subset\!\supset (\mathbf{f}\cdot\mathbf{\hat{n}})\,{\rm d}S = \Sigma \,\!

Equivalence between elementary and integral form

Starting from

 \iint_{S_1} \bold{f}_1 \cdot {\rm d}\bold{S}_1 = \iint_{S_2} \bold{f}_2 \cdot {\rm d}\bold{S}_2

the surfaces are equal (since there is only one closed surface), so S1 = S2 = S and we can write:

 \iint_{S} \bold{f}_1 \cdot {\rm d}\bold{S} = \iint_{S} \bold{f}_2 \cdot {\rm d}\bold{S}

The left hand side is the flow rate of quantity q occurring inside the closed surface S. This must be equal to

 \iint_{S} \bold{f}_1 \cdot {\rm d}\bold{S} = \Sigma - \frac{{\rm d}q}{{\rm d}t}

since some is produced by sources, hence the positive term Σ, but some is also leaking out by passing through the surface, implied by the negative term -dq/dt. Similarly the right hand side is the amount of flux passing through the surface and out of it, so

 \iint_S \bold{f}_2 \cdot {\rm d}\bold{S} = \iint\limits_{S}\!\!\!\!\!\!\!\!\!\!\!\subset\!\supset \bold{f} \cdot {\rm d}\bold{S}

Equating these obtains the integral form:

 \Sigma - \frac{{\rm d}q}{{\rm d}t} = \iint\limits_{S}\!\!\!\!\!\!\!\!\!\!\!\subset\!\supset \bold{f} \cdot {\rm d}\bold{S}
 \Sigma = \iint\limits_{S}\!\!\!\!\!\!\!\!\!\!\!\subset\!\supset \bold{f} \cdot {\rm d}\bold{S} %2B \frac{{\rm d}q}{{\rm d}t}

Electromagnetism

3-currents

In electromagnetic theory, the continuity equation can either be regarded as an empirical law expressing (local) charge conservation, or can be derived as a consequence of two of Maxwell's equations. It states that the divergence of the current density J (in amperes per square meter) is equal to the negative rate of change of the charge density ρ (in coulombs per cubic metre),

 \nabla \cdot \mathbf{J} = - {\partial \rho \over \partial t}

Maxwell's equations are a quick way to derive the continuity of charge.

4-currents

Conservation of a current (not necessarily an electromagnetic current) is expressed compactly as the Lorentz invariant divergence of a four-current:

J^\mu = \left(c \rho, \mathbf{j} \right)

where

since

\partial_\mu J^\mu = \frac{\partial \rho}{\partial t} %2B \nabla \cdot \mathbf{j}

then

\partial_\mu J^\mu = 0

which implies that the current is conserved:

\frac{\partial \rho}{\partial t} %2B \nabla \cdot \mathbf{j} = 0.

Interpretation

Current is the movement of charge. The continuity equation says that if charge is moving out of a differential volume (i.e. divergence of current density is positive) then the amount of charge within that volume is going to decrease, so the rate of change of charge density is negative. Therefore the continuity equation amounts to a conservation of charge.

Fluid dynamics

In fluid dynamics, the continuity equation states that, in any steady state process, the rate at which mass enters a system is equal to the rate at which mass leaves the system.[1][2] In fluid dynamics, the continuity equation is analogous to Kirchhoff's current law in electric circuits.

The differential form of the continuity equation is:[1]

 {\partial \rho \over \partial t} %2B \nabla \cdot (\rho \mathbf{u}) = 0

where

If ρ is a constant, as in the case of incompressible flow, the mass continuity equation simplifies to a volume continuity equation:[1]

\nabla \cdot \mathbf{u} = 0,

which means that the divergence of velocity field is zero everywhere. Physically, this is equivalent to saying that the local volume dilation rate is zero.

Further, the Navier-Stokes equations form a vector continuity equation describing the conservation of linear momentum.

Energy

By conservation of energy, which can only be transferred and not created or destroyed leads to a continuity equation, an alternative mathematical statement of energy conservation to the thermodynamic laws.

Letting

the continuity equation is:

 \nabla \cdot \mathbf{q} %2B \frac{ \partial u}{\partial t} = 0

Quantum mechanics

In quantum mechanics, the conservation of probability also yields a continuity equation. The terms in the equation require these definitions, and are slightly less obvious than the other forms of volume densities, currents, current densities etc., so they are outlined here:

 \rho = \Psi^{*} \Psi \,\!
P=P_{\mathbf{r} \in V}(t) = \int_V \Psi^{*} \Psi \mathrm{d} V = \int_V |\Psi|^2 \mathrm{d} V \,
 \mathbf{j} = \frac{\hbar}{2mi} \left [ \Psi^{*} \left ( \nabla \Psi \right ) - \Psi \left ( \nabla \Psi^{*} \right ) \right ].

With these definitions the continuity equation reads:

 \nabla \cdot \mathbf{j} %2B \frac{\partial \rho}{\partial t} = 0 \rightleftharpoons \nabla \cdot \mathbf{j} %2B \frac{\partial |\Psi|^2}{\partial t} = 0.

Either form is usually quoted. Intuitively; the above quantities indicate this represents the flow of probability. The chance of finding the particle at some r t flows like a fluid, the particle itself does not flow deterministically in the same vector field.

See also

References

  1. ^ a b c Pedlosky, Joseph (1987). Geophysical fluid dynamics. Springer. pp. 10–13. ISBN 9780387963877. 
  2. ^ Clancy, L.J.(1975), Aerodynamics, Section 3.3, Pitman Publishing Limited, London
  3. ^ Quantum Mechanics Demystified, D. McMahon, Mc Graw Hill (USA), 2006, ISBN(10) 0 07 145546 9

Further reading